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Abstract 

The emission from neutral hydrogen (HI) clouds in the post-reionization era < 6), too faint to 
be individually detected, is present as a diffuse background in all low frequency radio observations 
below 1420 MHz. The angular and frequency fluctuations of this radiation (~ 1 mK) is an important 
future probe of the large scale structures in the Universe. We show that such observations are a very 
effective probe of the background cosmological model and the perturbed Universe. In our study 
we focus on the possibility of determining the redshift space distortion parameter (5, coordinate 
distance r^, and its derivative with redshift, r'^. Using reasonable estimates for the observational 
uncertainties and configurations representative of the ongoing and upcoming radio interferometers, 
we predict parameter estimation at a precision comparable with supernova la observations and 
galaxy redshift surveys, across a wide range in redshift that is only partially accessed by other 
probes. Future HI observations of the post-reionization era present a new technique, complementing 
several existing one, to probe the expansion history and to elucidate the nature of the dark energy. 
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I. INTRODUCTION 



Determining the expansion history of our Universe and parameterizing the constituents 
of the Universe at a high level of precision, are currently some of the most important goals 
in cosmology. While high-redshift {z < 2) supernova la observations (e.g. [1, 2]) and galaxy 
surveys < 1 ) (e.g. [3]) probe the local universe; and CMBR observations (e.g. [4, 5]) 
probe the recombination era {z ~ 1000), the expansion history is largely unconstrained 
across the vast intervening redshift range. Observations of redshifted 21cm radiation from 
neutral hydrogen (HI) hold the potential of probing the universe over a large redshift range 
(20 > z > 0): from the dark ages to to the present epoch (eg. [6, 7]). Such observations 
can possibly be realized at several redshifts, using the currently functioning GMRT [34]. 
Several new telescopes are currently being built with such observations in mind (eg. MWA 
[35] & LOFAR [36]). Such observations will map out the large-scale HI distribution at high 
redshifts. It has recently been proposed [8, 9] that Baryon Acoustic Oscillations (BAO) 
in the redshifted 21cm signal from the post-reionization era {z < 6) is a very sensitive 
probe of the dark energy. The BAO is a relatively small (~ 10 — 15 per cent) feature that 
sits on the HI large-scale structure (LSS) power spectrum. In this paper we investigate 
the possibility of probing the expansion history in the post-reionization era using the HI 
LSS power spectrum without reference to the BAO. Unless otherwise stated we use the 
parameters {flmo, ^ao, ^bh^, h, Us, erg) = (0.3, 0.7, 0.024, 0.7, 1.0, 1.0) referred to as the LCDM 
model in our analysis. 

At redshifts z < 6, the bulk of the neutral gas is in clouds that have HI column densities 
in excess of 2 x 10^° atoms/cm^ [10-12]. These high column density clouds are observed 
as damped Lyman-a absorption lines seen in quasar spectra. These observations indicate 
that the ratio of the density Pgas(-2) of neutral gas to the present critical density Pcrit? of 
the universe has a nearly constant value Pgas(-2)/Pcrit ~ 10~^, over a large redshift range 
< z < 3.5. This implies that the mean neutral fraction of the hydrogen gas is xm — 
50 Qgas/i^(0.02/Q6/i^) = 2.45 X 10~^, which we adopt for the entire redshift range z < 6. 
The redshifted 21cm radiation from the HI in this redshift range will be seen in emission. 
The emission from individual clouds (< 10 fiJj) is too weak to be detected with existing 
instruments unless the image is significantly magnified by gravitational lensing [13]. The 
collective emission from the undetected clouds appears as a very faint background in all 



2 



radio observations at frequencies below 1420 MHz. The fluctuations in this background with 
angle and frequency is a direct probe of the HI distribution at the redshift z where the 
radiation originated. It is possible to probe the HI power spectrum at high redshifts by 
quantifying the the fluctuations in this radiation ([14, 15]). 



II. FORMULATION 

The Multi-frequency Angular Power Spectrum ( MAPS) Gi{lS.v) [16] quantifles the statis- 
tics of the HI signal as a joint function of the angular multipole i and the frequency separa- 
tion 1S.V. We deflne the angular power spectrum Cn — C^(0) and the frequency decorrelation 
function 

M-)^^. (1) 

to separately characterize the angular and the Ai/ dependence respectively. The latter 
quantifles whether the HI signal at two different frequencies v and v -|- Ai/ is correlated 
K^(Ai/) ~ 1 or uncorrelated K^(Ai/) ~ . The function C^(Ai/) can be estimated directly 
from observations without reference to a cosmological model (eg. [17]). However, it is 
necessary to assume a background cosmological model in order to interpret CelAu) in terms 
of the three dimensional LSS HI power spectrum. On the large scales of interest here, it is 
reasonable to assume that HI traces the dark matter with a possible linear bias b, whereby 
the three dimensional HI power spectrum is b'^P{k), where P{k) is the dark matter power 
spectrum at the redshift where HI signal originated. We have [16] 

Ci{Au) = — / d^ii cos(^|| rl Au) Pni{k) , (2) 
^'"i^ Jo 

where the three dimensional wavevector k has been decomposed into components k\\ and 
I/Ti,, along the line of sight and in the plane of the sky respectively. The comoving distance 
Tj, is the distance at which the HI radiation originated. Note that (1 -|- z)~^ = d^iz) is the 
angular diameter distance and — drv/dv. The temperature occurring in eq. (2) is given 



^0.02 J \ h J H{z) 

and PHi(k) is the three dimensional power spectrum of the "21 cm radiation efficiency in 
redshift space" , which in this situation is given by 

Pm(k)=4i&Ml+/5/^')'^W- (4) 
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FIG. 1: Here we plot C;(0) at redshifts z = {1.5, 3.0, 4.5}. The signal decreases monotonically with 
increasing redshift, so the lowest plot is for the highest redshift. We assume the bias to be 6 = 1 
throughout. 

The term (1 + arises due to HI peculiar velocities ([14, 18]), which we assume to be 

determined by the dark matter. This is the familiar redshift space distortion seen in galaxy 
redshift surveys, where n = k\\/k. and P = f{z)/b is the linear distortion parameter, which 
is the ratio of f{z) that quantifies the growth rate of linear perturbations, and b the linear 
bias. 

III. RESULTS AND CONCLUSIONS 

The expected signal C/(Ai/) from a few representative redshifts, calculated for the 
LCDM model, is plotted in Figure 1, and in Figure 2 we have plotted the frequency 
decorrclation function fi;^(Az/) as a function of Au, for a fixed redshift z = 3.0 and for 
£ — 100, 1000 & 10000. The HI signal is smaller than ~ 1 mK, and it decreases with in- 
creasing I. The shape or £ dependence is decided by the shape of P{k) at all comoving 
wave-numbers k > £/r,^. The signal at two different frequencies u and u + Av decorrelates 
rapidly with increasing Au and K^(Ai/) < 0.1 at Au > 5 MHz. The decorrelation occurs 
at a smaller Au for the larger multipoles (Figure 2). While the HI signal at a frequency 
separation Au > 5 MHz is expected to be uncorrelated, the foregrounds are expected to be 
highly correlated even at frequency separations larger than this (eg. [19]). This should in 
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FIG. 2: Here we plot the frequency decorrelation function kh^/S.v) as a function of Ai^, for a fixed 
redshift z = 3.0 and i = {100, 1000, 10000}. The signal declines more sharply for higher value of £. 



principle allow the HI signal to be separated from the foregrounds, which are a few orders 
of magnitude larger (eg. [20, 21]). 

It is clear from eq. (2) that Ce{Ai/) depends on the background cosmological model 
through the parameters r,^, r'^). Assuming that the dark matter power spectrum P{k) is 
known a priori, observations of C^(Ai/) can be used to determine the values of these three 
parameters. It is convenient to replace with the dimensionless parameter [22] 

Figure 3 shows the variation of the three parameters ri,,p) across the redshift range z < 6 
for the LCDM model. 

We separately consider parameter estimation using Q and Ki{Ai/). The former does not 
depend on p. The amphtude A — [TxmbY / t^t^ of Q is uncertain, and we consider the 
joint estimation of three parameters {A, /3, r^) from observations of C^. The value of K(,{/S.v) 
is insensitive to the amplitude A, leaving three parameters {P,r^,p) that can be jointly 
estimated from this. We use the Fisher matrix (e.g. [23]) to determine the accuracy at 
which these parameters can be estimated. 

Parameter estimation depends on two distinct aspects of the observing instrument. The 
first is the £ range ie. imin, ^max, and the sampling interval Ai, which corresponds to the 
smallest i spacing at which we have independent estimates of Q( Ai/) . This is determined by 
the instrument's field of view, and is inversely related to it. The second is the observational 




FIG. 3: Here we plot the parameters (/?, r,p) as a function of redshift z for the concordance LCDM 
model. The parameter r = r,^/(6000 Mpc). 

uncertainty in Cei^Au). This is a sum, in quadrature, of the instrumental noise and the cosmic 
variance. The cosmic variance contribution 6Ce/Ce = ^2/{{2i + 1) i Ai) (f is the fraction of 
sky observed) is further reduced because the large frequency bandwidth Ai^b provides several 
independent estimates of Ce. We assume that 5Ci is reduced by a factor we ^/ Aub/{1 MHz) 
because of this. The instrumental uncertainties were estimated using relations [17] between 
SCi and the noise in the individual visibilities measured in radio- interferometric observations. 
For this we assume that the baselines in the radio-interferometric array have a uniform u-v 
coverage. 

We consider three different instrumental configurations for parameter estimation. 

A. The currently functional GMRT has too few antennas for cosmological parameter 
estimation. We consider an enhanced version of the GMRT with a substantially larger 
number of antennas {N — 120) , each identical to those of the existing GMRT. The 
antennas have a relatively small field of view (^fwhm ~ 0.8° at 610 MHz) and the array 
has relatively large baselines spanning = 500 to £jnax — 10, 000 with A£ — 100. 

B. The upcoming MWA will have a large number of small sized antennas. The antennas 
have a relatively large field of view (^fwhm ~ 5° at 610 MHz), and the array is expected 
to be quite compact spanning = 100 to ijnax — 2000 with Ai — 20. The first 
version of this array is expected to have N — 500 antennas which is what we consider. 
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FIG. 4: Expected one-sigma fractional errors for parameter estimation at different redshifts for the 
LCDM model. The curves in each panel correspond, from top to bottom, to the cases A, B, and 
C, respectively. 

C. This is a future, upgraded version of the MWA which is expected to have N — 5000 
antennas. 

For each of these configurations, we assume that 16 simultaneous primary beams can be 
observed. We present results for 2 years of observation for A and B, and 1000 hours for 
C. Throughout we assume frequency channels 0.05 MHz wide, a bandwidth Aub — 32 MHz, 
and that a single field is observed for the entire duration. For parameter estimation we use: 
6Ki{Au) = V26Ce/Ce. 

We find that observations of impose very poor constraints on the parameters /? and 
r^, and we do not show these here. The accuracy is considerable higher for «;^(Ai/), which 
captures the three dimensional clustering of the HI as compared to Q, which quantifies 
only the angular dependence. Figure 4 shows the predicted estimates for the parameters 
(3, and p at various redshifts. Further, we find that a compact, wide-field array (B,C) is 
considerably more sensitive to these parameter as compared to case A. 

Considering the three parameters individually: 

Redshift-space distortion parameter: /3 . This has traditionally been measured from 
galaxy redshift surveys [24-27], with uncertainties in the range 0.1 < A/3//3 < 0.2. These 
observations have, till date, been restricted to z < 1. Future galaxy surveys are expected to 
achieve higher redshifts and smaller uncertainties. Galaxy surveys have the drawback that 
at very high redshifts they probe only the most luminous objects, which are expected to be 
highly biased. HI observations do not have this limitation and could provide high precision 
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FIG. 5: Expected one-sigma confidence regions for the parameters QmO and U^q, based on estimated 
errors for observations of p, corresponding to Figure 4, at z = 3. 

(A/3//3 < 0.1) estimates over a large redshift range. 

Coordinate distance, r,^: The most direct measurement of the coordinate distance 
comes from supernova type la observations for 2; < 2. Current Sn la observations give 
^Ty/ri, ~ 0.07 [28] for a single supernova. The statistical error in the coordinate distance 
can be further reduced by observing a large number of supcrnovac in a small redshift bin; 
thus the fundamental limitation of this technique is due to unknown systematics in the 
supernovae themselves, since it is certainly possible that supernovae at high redshift are 
different. Figure 4 shows that the HI method might have the potential to enable a precise 
measurement of the coordinate distance up to much larger redshifts. Furthermore, such a 
complimentary probe will also help in ascertaining systematics in the supernova probe. 

Derivative of coordinate distance, p: This quantifies the Alcock-Paczynski (AP) 
effect [29], which is well accepted as a means to study the expansion history at high 
though such observations have not been possible till date. Observations of redshifted 21cm 
radiation hold the potential of measuring the AP effect [22, 30, 31]. The parameter p is not 
affected by the overall amphtude A and the bias 6, and is a sensitive probe of the spatial 
curvature (Figure 3). Our estimates indicate that it will be possible to measure p with an 
accuracy /\p/p ~ 0.03 over a large z range. 

The parameters {j3,ry,p) chosen for our analysis occur naturally when we interpret 
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C^(Ai/) in terms of the three dimensional dark matter power spectrum P{k). Further, 
these parameters are very general in that they do not refer to any specific model for either 
the dark energy or the dark matter, and are valid even in models with alternate theories of 
gravity (eg. [32, 33]). In fact, observations of these three parameters at different redshifts 
can in principle be used to distinguish between these possibilities. 

For the purpose of this paper, we illustrate the cosmological parameter estimation by 
considering the simplest LCDM model, with two unknown parameters fi^o and Qfeo, and 
^Ao — ^ ~ ^mo — ^ito- In Figure 5 we plot the 1-a confidence interval for the estimation of 
ftrno and Qfeo, using a single measurement of p alone, ie. only one of the three parameters 
measured at a single redshift z = 3. Note that p is insensitive to Hq and hence it is not 
considered as an additional parameter here. It is possible to combine measurements at 
different z to improve the constraints on cosmological parameters. We shall undertake a 
detailed analysis for quantifying the precision that can be achieved by combining different 
data sets (CMBR, galaxy surveys) for a more comphcated dark energy model in a future 
work. 

In conclusion, HI observations of the post-reionization era can, in principle, determine the 
expansion history at a high level of precision and thereby constrain cosmological models. 
Neither the upcoming initial version of the MWA which is planned to have 500 antenna 
elements nor any conceivable upgradation of the existing GMRT will be in a position to 
carry out such observations, the observation time needed being too large. We find that 
an enhanced version of the MWA, which is planned to have 5000 antenna elements, would 
be in a position to meaningfully constrain cosmological models. By combining different 
probes, we expect to achieve an unprecedented precision in the determination of cosmological 
parameters. This will be a step towards pinning down the precise nature of dark energy in 
the universe. 
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